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Abstract 

In this work we investigate the general properties and the ground state of an asymmetrical dilute 
gas of cold fermionic atoms, formed by two particle species having different densities. We have 
shown in a recent paper, that a mixed phase composed of normal and superfluid components is the 
energetically favored ground state of such a cold fermionic system. Here we extend the analysis 
and verify that in fact, the mixed phase is the preferred ground state of an asymmetrical superfluid 
in various situations. We predict that the mixed phase can serve as a way of detecting superfluidity 
and estimating the magnitude of the gap parameter in asymmetrical fermionic systems. 
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I. INTRODUCTION 



The recent experimental advances in cold atomic traps have induced a great amount 
of interest in fields from condensed matter to particle physics, including approaches and 
prospects from the theoretical point of view. The advent of new techniques in dealing 
with cold Fermi atoms in traps in the last few years [3, IE IE HI IE IE 0] have revived the 
interest in some fundamental questions in many-body quantum physics. The ability of 
preparing a trap with a strongly interacting ultra-cold gas of fermionic atoms 0] constitutes 
an example of such an experimental development. In this scenario one may found the 
necessary conditions for pair formation between species whose Fermi surfaces do not match. 
An immediate question, which arises in this context, is what is the ground state of a fermion 
system composed by two components having different masses and densities (or chemical 
potentials)? 

In this paper we study the possibility and the prerequisites for pair formation between two 
species of particles with unequal number densities or chemical potentials. These different 
species could be two fermionic atoms ( 40 K or 6 Li) 0, 0, 0] or hyperfine states of the same 
atom 0]. The nucleus of neutron stars could also have the basic ingredients for asymmetric 
pairing, since quarks with different flavors would have different densities |E EE n] • 

The pairing between particles with opposite spin and equal and opposite momenta near 
their (common) Fermi surface, is described by the standard Bardeen-Cooper-Shrieffer (BCS) 
T2l | theory of superconductivity. A model exhibiting asymmetry between the (fixed) chem- 



ical potentials of two particle species was first studied by Sarma a long time ago at zero 
and finite temperature ^E]. The Sarma phase was proposed recently as a candidate for 
the ground state of a system of two species of cold fermionic atoms 0] . This state of 
matter, named "internal gap" phase, would be composed by a (homogeneous) mixture of 
normal and superfluid phases. 

It has been showed in a recent paper ^E] that the ground state of such an asymmetrical 
fermionic system is: I. either the normal or the BCS states, depending on the asymmetry 
between the Fermi surfaces of the two species, for the situation of fixed chemical potentials, 
II. a (inhomogeneous) mixed phase formed by normal and superfluid components, for the 
case of fixed number of particles and finally III. a mixed phase having as limits the BCS 
state for m a = and normal phase with only b particles when mj 3> m a , for fixed overall 
number density. 

It has been shown [IE] that the coupling constant in the internal gap (Sarma phase) has 
a lower value, g m ini instead of being arbitrarily small as in the BCS phase. We show in 
this paper that the coupling constant supported by the Sarma phase also has a maximum 
value, g ma x- Thus, the Sarma phase can exist only in the window g m in < g < g-max- As a 
consequence, the Sarma phase is not favorable for weak enough attractive interactions. For 
the sake of comparison, we show that the approximate number of particles available to form 
pairs in the Sarma phase is always less than the one in the BCS phase. As a result, the 
Sarma phase may have a disadvantage in the condensation energy, since the energy of the 
system is lowered with pairing. Even for the case of fixed chemical potentials, we found a 
particular combination for the Fermi surfaces P F and P F of the two species, for which there 
exist a mixture of normal and BCS phases. In the case of fixed number of particles, we 
found that the state of lower energy satisfies the conditions for equilibrium between the two 
components of the mixed phase. 

We notice that another possibility for the ground state in asymmetric matter is the Larkin, 
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Ovchinnikov, Fulde and Ferrel (LOFF)-phase The gap parameter in the LOFF phase 
also breaks translational invariance, as the mixed phase does. However, since the LOFF 
phase can exist only in a very narrow window of asymmetry for the chemical potentials we 
do not consider this possibility in our discussion. 

The counterpart for ordering in quantum chromodynamics (QCD) (for the case where 
the quarks have approximately the same Fermi surface) at high baryon density is the color 
superconductivity 0, EH EH E3 • For high density asymmetric quark matter the analog of 
LOFF state leads to crystalline color superconductivity I n Q CD with two light flavors, 
the ground state is the two-flavor color superconductor (2SC) |2Cj. 

Still in quark matter, it has been argued that the gapless color superconductivity (anal- 
ogous to the Sarma phase) can reach stability, under the conditions of color and electrical 
charge neutrality, in the two j^, 2^1 and three |23[ flavor cases. 

The paper is organized as follows. In Section II we introduce the model Hamiltonian and 
derive the free energy and gap equation by variational and mean field approximations. In 
Section III we search for the fundamental state of the asymmetrical fermion superfluid in 
three situations: fixed chemical potentials, fixed particle numbers and fixed overall density. 
We discuss the results and conclude in Section IV. 



II. THE MODEL, BASIC DEFINITIONS AND THE GAP EQUATIONS 

To begin with, let us consider a nonrelativistic dilute (i.e., the particles interact through 
a short-range attractive interaction) cold fermion system, described by the following Hamil- 
tonian 

% = H - ^ H<* n a = ^2 £ k a k a k + £ k b k b k +9^2 a k' b -k' b -k a k, (1) 
k,a k k,k' 

where a k , are the creation and annihilation operators for the a particles (and the same for 
the b particles) and e k are their dispersion relation, defined by e k = — fi a , with = ^— 
and /x a being the chemical potential of the a-specie, a = a,b. To reflect an attractive 
(s-wave) interaction between particles a and b we take g < 0. We shall derive below the 
state encountered by an asymmetric fermion superfluid at zero temperature, i.e., when its 
energy is a minimum. This investigation shall be done in two scenarios representing different 
physical situations, namely fixed chemical potentials and fixed number of particles. 



A. Solution by a variational procedure 



The BCS ground state, which describes a superposition of empty and occupied (paired) 
states, is given by 



\bcs) = n [ 



u k + Vkal^_ k 



|0> 



(2) 



where the arbitrary complex (a priori) coefficients Uf. and Vk are to be determined by a 
variational calculation. They are subjected to the normalization, |u^| 2 + |t^| 2 = 1, and spin 
singlet, Uk = U-k, Vk = f-fe, conditions. These constraint are satisfied if one sets u k = sin^ 
and Vk = cos 9k III • 
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The free-energy (or thermodynamical potential) for this system is obtained as usual by 
taking the expectation value of !K 

F =< BCS|:K|BCS >= ^4cos 2 fl fc + !^Bin20 fc sin20 fc /. ( 3 ) 

k,a k,k' 

When the minimum of F with respect to 6 k is required (M- = 0), we get tan 26 % = g ^ k ' S1 ^ 26,fc/ ; 

k 

where e~£ = k 2 k = 2 ( 2M — ^ ) ' w ^h M = ^ a ^ b being the reduced mass and \x = fi a + /•*&• 
Defining the gap parameter as 

A = -§£ S in20 fc ,, (4) 
1 v 

it follows that tan 20*. = — This last equation combined with the normalization condition 

e k 

gives 

ut * = l( i+e 4) =i -*> (5) 



where E = y e^ 2 + A 2 . Substituting these quantities back in Eq. @ we obtain the final 
expression for the free-energy 



(2tt) 3 k V £ 



together with the gap equation 
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d 3 k 



1 ~2j r^— t - < 7 > 



Following Ref.|25|], the gap equation can be related to the two-body scattering length a, in 
dimensional regularization technique, and we have 

M If d 3 k 1 

(8) 



2na 2 J (2vr)3 



The change of variable = z allows us to write J — > A ^j F Jo° Z2 dz, where fc? 
a/2~M/j, so that 



7T 



fci^a \ [i ) 

where 



-A/2 ( — ) , O) 



^ 7T 



= <** „ = . „ = --fsza + *r /2 Pe ), (io) 



^{z-l) 2 + x 2 sm/3vr v V v 7 ! 



X 2 
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with P/3 being the Legendre function. In the weak coupling limit, A/u << 1, we expand 
the function above and obtain, to leading order, 



. 4u 

A ~ f>2kp\a\ 

e 2 

The free-energy, Eq. (JBJ), can be obtained in the same way by the use of Eq. (jl(Jj) : 
A 2 Mfi 2 k F 



F = — 

9 



4vr 2 



(11) 



(12) 



U5 
F 



30tt 2 M 2tt 2 
For the number densities we have 



Mk F A 9 
A 2 . 



dF 



ii, 



dF dn 



1.3 



M 2 A 2 



<9F 
9// 



(13) 



+ 



67T 2 2lX 2 kp 



71 



\a\kr 



n. 



Then the particles a and b have the same density in the BCS phase, as it should be. The 
energy of the superfluid state is 



E = F + ^u a n a = F + 



(14) 



b 5 
n p 



+ 



klA 2 



20n 2 M 8vr 2 u 



7T 



We can invert Eq. (fT3|) expressing u = f{n) so we can write the energy as a function of the 
number density 



E(n) 



(Qn 2 n 



2„\5/3 



M 



20tt 2 M 2tt 2 
(67r 2 n) 5 / 3 
20tt 2 M 



(15) 



]_ _ 4Q e |a|(67r2„)l/3 



E BCS (n). 



Taking the difference in energy between the superconducting and the normal state we find 
the condensation energy 



E BCS (n) -E N (n) 



M 



2tt' 



(Gt^^A 2 



(16) 



B. Solution by the (Bogoliubov) canonical transformation 

Since we are considering a dilute system, we can make the mean-field approximation 
24l 26[. Defining A = —g J2k < ^-k a k > that is chosen as being real, we get 



A 2 

ft = Yl e >i ak + e * 6 * 6 * - A ( a i b ~k + - — • 

1. y 



(17) 
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QPE 




FIG. 1: Dispersion relation for the quasi-particles a and j3 showing a region where £^ is negative 
for m b = 7m a , Pp = 1.45 PJ;, and A = 0.29A , obtained by Eq. Solid curve correspond to 

£^ and dashed curve correspond to 



Since now U£ is quadratic, it can be diagonalized 



w = E £ * a ^ + £ ^ + £ * - £ " - V' (18) 



where £^ ,/3 are the quasi-particle excitations (QPE) 

= ±e~ k +E = ±e- + y/ef + A*, (19) 

where e k = k 2 k . The (linear) transformation used in the diagonalization of "K is defined 

as 

' Oi k \ (u k -v k \ ( a k \ 



with Uk and given by Eq. (JHJ). Depending on the relative magnitudes of P F , P F , nib 
and m , either £^ have two roots, where Pp = (2m a /i Q )a are the particles Fermi surfaces. 
We chose P F > Pp and m;, > m a such that only £jj! crosses zero. The roots of £f can be 
determined from the equation 

e a A = -A 2 . (21) 

They are define as k\ and k 2 , 

^i 2 2 = + if) T \[{P b F ~ Pff ~ 16m a m 6 A 2 ] 1/2 - (22) 

Fig. ^ shows that, for some values of A, £^ may be negative for momenta k\ < k < k 2 . A 
simple inspection in fc 12 shows that A has a critical (maximum) value, above which /c 12 are 
not real, 

A. = J^L, (23) 

where 5p F = P£ 2 - Pf. 

The consequence of this is that we can have two physically distinct situations: 
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1. Symmetrical superfluid. 



The first situation is when A > A c . In this case £^ are both positive and we have the 
standart BCS with Eq. (flHjl reduced to 

K ^ ^bcs = J2e+-E-— J (24) 
k 9 

with the gap parameter given by Eq. (jllj) . The "average" Fermi surface for the asymmetrical 
BCS superfluid can be rewritten as 

k F = lpf + —^^8 PF . (25) 
V m a + m b 

Although we could write the gap with kp depending on 5pp by the equation above, it is a 
not increasing function of 5pp (with Pp, m a and m& kept fixed). As it will become clear 
below, the BCS free-energy minimum does not depend on the asymmetry between the Fermi 
surfaces of the two species until a point where there is a first order phase transition to the 
normal phase. 



2. Asymmetrical superfluid. 

The second situation, that we call "Sarma phase" (first pointed out in Ref. (13)) happens 
when A < A c . In this case £f < between fc 12 . The free energy in the Sarma phase 
is obtained when we find a state which minimizes the diagonalized Hamiltonian, Eq. 
f!18|) . The smallest energy is reached when the modes with negative ££ are filled and the 
remaining modes are left empty. More precisely, the ground state satisfies 

a k ,f3 k \V) = if £f>0, 

441*) = if £^<0. (26) 

This state can be written in terms of the original a\ and b k operators and the vacuum state 
|0) as 

m= n h+^4^ fc ]n^i°)- ( 2? ) 

k<ki k\ 
k>k 2 

The state above corresponds to having BCS pairing in the modes k where £^ ,/3 > and a 
state filled with particles b (a) in the modes where £^ < (£^ < 0). The free-energy of the 
Sarma phase turns out to be 

^^ s =L^-v+C <§rf-T- (28) 



k>k 

with = (remembering that the partial derivative also hits k\ and k2 in the limits of 
the integrals) we obtain the gap equation 

g I* d 3 k 1 

~2 Jk^ (2vr) 3 n^~r:' ^ ' 

J k>k 2 y > del + A 
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Again, the standard BCS result can be recovered for A > A c , in which case the equation 
above reduces to Eq. (JJJ) and the free energy goes to 




The equation above can be written in the form of Eq. (jUJ) using the gap equation (|7|). since 
Eq. is valid only at the minimum, due to the variational method we employed in the 
previous subsection. 



3. Solving the gap equation 



We show next that the solution for the gap equation in the Sarma phase can exist only 
for g min < g < g m ax- Since the integral in Eq. (j2Tlj) is dominated for momenta around kp, 
let us restrict ourselves to the region (k F — A) < \k\ < k\ and k 2 < \k\ < (kp + A), where 
A is an ultraviolet cutoff (used only for the sake of this proof, none of our results depend 
on it) and suppose, for this analysis, that P F > P F [3, 



grain is found when A has its 



maximum value namely, A r 



iV(0)ln 



2 A kp (rn a + m b ) 
5p F ^m a m b 




2 A kp (m a + m b 



5p F ^m a m b 



-i 



(31) 



where N(0) = J -^^S(e^ ) = ^f- is the density of states. The maximum value of g is 
obtained when A = 0, 



7V(0) In 



2 X kp {m a + m b ) 
5p F y/m a m b 



A consequence of this is that for 2 , A kp ( m "±!^) » 1 (or ^ 
g m i n and g max is constant and independent of the cutoff A 

= -ln2iV(0). 



1^ 



\9r, 



(32) 

oo) the window between 
(33) 



Let us denote the gap of the standard BSC, which is given by Eq. (jlljl . by Ao- We can find 
the gap in the Sarma phase through the identity 



M 
27r|a| 



d 3 k 



(27T) 



A 2 



d 3 k 



fc<fci (27r) 



k>k 2 



(34) 



A s 



For small values of the gaps (A , A << fi a , fi b ) the integrals can be approximated and it is 
found that 

A 2 eKh) 



A 2 



This equation for A can be expressed as 

A 2 _ Sp F - [Sp 2 F - lGm.mfeA 2 ] 1 / 2 
A 2 ~ 5pp + [5p 2 F - 16m a m 6 A 2 ] 1 /2 ' 



(35) 



(36) 



S 



It is easy to see that when A = A r we have one solution, Ag = A c = Ao, which is the BCS 
case. The other solution for A is k| 27J 



A*^ WAo(-^L-Ao). 



2y/m a m b 

This equation for A generalizes the Sarma results The conclusion is that 

2^m a m b A < \5p F \ < 4 v /m a m 6 A (38) 
<A 5 < A . 

The number densities in the Sarma phase are given by n a = — ^p— , then we get 



n b = n a + 



k>k2 v ^" ; \ \l e%* + A 2 
fc2 d 3 k 



ki (2vr) 3 

It shows that outside the region k\ < k < k 2 the particles a and b are paired via BCS having 
number density n a , whereas inside the region of momenta k\ < k < k2 there are only b 
particles. We define the particle occupation numbers for the species a and b respectively, as 




u a = J I 1 £k I , for k < ki and k > k 2 , (40) 



ef + A 2 



u b = - [ 1 . k = I , for k < k\ and k > k 2 , and 1 for k\ < k < k 2 . 



2 



et + A 2 



In Fig. @ we show the particle occupation numbers for the species a and b. We also plot 
the hierarchy of momenta for those values of the parameters we have used in Fig. As 
we can see, the number of particles of both species having momenta k > k 2 is negligible. 
This permits us to say that the approximate number density of particles available to form 
pairs in the Sarma phase is 

This fact can have consequences in the condensation energy of the Sarma phase since, as we 
have seen in Eq. (|16|). pairing lowers the energy of the system. 



III. SEARCHING FOR THE LOWEST ENERGY 
A. Fixed chemical potentials 

The situation where the chemical potentials are kept fixed can find place for instance in 
a gas of fermionic atoms connected to reservoirs of species a and b so the number densities 
are allowed to change in the system. 
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FIG. 2: Figures on top: On the left it is shown the particle occupation number u a for the specie a 
in the Sarma phase as compared with the BCS state (dashed curve). On right there is a amplified 
plot of the specie a occupation number for k > ki . Figures on bottom: On left it is shown the 
particle occupation number u b for the specie b in the Sarma Phase also compared with the BCS 
state (dashed curve). Again we show on right a amplified plot of the specie b occupation number 
for k > &2. We also show the hierarchy between the momenta involved, P F < k\ < kp < k2 < Pp. 



In Fig. El we show the thermodynamic potential F s as a function of A for different values 
of Pp and Pp, keeping the combination kp/2M = Pp 2 /2m a + P F /2m& fixed, computed 
from a numerical evaluation of Eq. (|28|). When As is larger than A c = ^=±= , k\ 2 are not 

real, the quasi-particle excitations £^ ,/3 are always positive, and the free-energy is unchanged 
from the Pp = P F situation. However, for As < A c , k\ 2 are positive and the free-energy can 
be lowered by filling the states between k\ and k 2 with b particles. The A = A corresponds 
to the BCS phase. The normal state (A = 0), in particular, can have its free-energy lowered 
with increasing P F — P F , and after P ff >m ~ ( see discussion on Pf,dm below) becomes 
smaller than the previous minimum at A . Then, for fixed kp and increasing Pp — P F , there 
is a first-order phase transition between the sup erfluid and the normal state, as it has been 
found in various physical situations jH, 2^, 3cf . 

We can understand this results with very simple physical terms in the following way. 
Starting with \x a = we have the BCS state with its free-energy given by Eq. (fPfj) 

F(/i a = n h = n) = E(n) - (fi a + fi b )n = E(n) - fi To tai n, (42) 

k 2 

where /i Totai = 2/i = ^§ is kept fixed. 

In order to compare the energies of the BCS state with another state having \i h > ji a 
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FIG. 3: thermodynamic potential F s for different values of P F and P F (kp fixed). The top curve 
corresponds to Pp = P F and the lower curves correspond to increasing values of \P F — P F \. 

(P F > Pp), we let the system absorb one b particle from the 6-particles reservoir 1 and 
eliminate one a particle to the a-particles reservoir (Again this will increase the energy of 
the system by at least Aq, by breaking a pair and promoting the remaining b particle to 
the elementary excitation levels). Then, the free-energy of the system having unequal Fermi 
surfaces turns out to be 

F(fi b > fi a ) = E(n) + 2A - fi a (n - 1) - fi b (n + 1) (43) 

= F(fJ, a = fl b ) + 2A - (fj, b - fl a ), 

where we have used that fi a + fi b = l^Totai- By absorbing one b particle and eliminating one 
a particle the system reduces its free-energy due to the (fib — fJ-a) term but also increases F 
due to the 2Ao term. Then, this is energetically favorable only if the difference in chemical 
potentials is large enough or the gap is small enough. In spit of the mismatch in the chemical 
potentials, the BCS state with equal number of a and b particles remains the ground state 
while 

A > 5fi, (44) 

where 5fi = B^-. 

As one can sees in Fig. El there exist a special combination of Pp and P F for which F s 
has double minima. Then, given a P F , we want to know the correspondent value of Pp that 
satisfies this requirement. We note that when A = — > k\ = P F and &2 = Pf whereas for 



Pp is 



p% 2 +p b2 



A — > h = k 2 = y F 2 F ■ Thus, the condition to find the relation between P F and 

F(A = 0) = F(A = A ) (45) 
d 3 k . , . ... [ p F d 3 k b f d 3 k . , A 2 

'4 - 141) + / 7S3s4 = / iw-^M - E )-—- 



,<P f(27r)3 V k I kl^J (27r) 3 * J (27r)3 V k 



1 This will increase the energy of the system by at least the value of Ao . The extra b particle added will be 
located at the level of elementary excitations. 
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FIG. 4: Range of momenta allowed for the excess of particles of specie b in the Sarma phase, where 
< A < A c . 



Since F is being evaluated at its minimum, we can substitute the r.h.s. of the equation 
above by Eq. JBJ) (whose solution is given by Eq. (O) and the l.h.s. is just the free-energy 
of the unpaired a and b particles 

1 

The solution of the equation above we cal PpDMi where DM stands for double minima. We 
point out here that this situation means coexistence of phases. Given these fixed values of 
Pp and Pp (or \i a and found above, the system could either be found in the normal, 
BCS or in a mixed phase, since those states would have the same minimal energy. When the 
free-energy F(P FDM ) has its minimum at Ao, after a increase in P F the gap jumps from Ao 
to 0, characterizing a first order phase transition from the superfluid to the normal phase, 
as we have discussed already. 

According to Eq. ()38|) there is (for some values of the Fermi surfaces Pp and P F ) an 
unstable state, that is the Sarma phase we have discussed in the previous section, in addition 
to the stable (or metastable) normal and BCS phases. This state corresponds to a maximum 
of F as a function of A situated between the BCS minimum A and the normal phase at 
A = 0. 

We can cast some light in the understanding of Figures El and |3] (and consequently learn 
what really happens in the situation of fixed chemical potentials in the Sarma phase (0 < 
A < Ac))- With that purpose we define a quantity D(A) = k 2 — k\ which is the region of 
momenta allowed for the (unpaired) b particles in the Sarma phase. For any curve having 
momenta Pp > Pp (and fixed in each curve) from the second curve of Fig. El until P ff> m 
the minimum of the free-energy of the system is achieved in the following way. When the 
(weak) attraction is turned on and increases, the system keeps releasing the "excess" of b 
particles until the point D(A C ) = (see Fig. HJ) where the BCS condition is reached and 

both species have the same density n a = n b = ^ + 27r2fc ^ ( 5 + J ■ For Pp > P F DM the 

minimum of F is obtained when the attractive interaction between the fermions decreases 
until it is turned off (D(A = 0) = Dm = Pp — Pp)- There is absorption of particles from 
the reservoirs so the asymmetry in the number densities gets its maximum and the system 

is in the normal state with number densities n a = -A- and rib = i^t- 



pa 5 



+ 



yb 5 



lc 5 

h,p 



Mk f 



30tt 2 M 2tt 2 



■A 2 



(46) 
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B. Fixed particle number 



The analysis regarding the ground state of asymmetrical fermion superfluids we had so 
far is valid for systems where the chemical potentials of the different species are kept fixed. 
We turn now to the cases in which the particle numbers are the fixed quantities. This is the 
relevant situation for trapped cold fermionic atoms. As we have seen, in the BCS state the 
particles have the same density, whereas in the normal ground state the number densities 
can be different. A mixed (inhomogeneous) phase has been proposed in Ref. 0] as being 
composed by islands of normal matter in a sea of BCS phase. According to Eq. (}3*9*j) the 
Sarma phase also supports different particle numbers, depending on the value of the gap. 
In order to find the ground state of a system with fixed n a and rib number densities it is, 
then, natural to compare the energies of the normal, Sarma and mixed phases. 



1. The mixed phase energy 

In the mixed phase, n a and rib particles in a box with volume V, are accommodated in 
such a way that in a fraction x of this volume the particles are free having densities n a} b, 
and in the rest of the volume there is pairing formation between the species a and b with 
number densities n BCS = nf cs = n. Then, the number densities in each component of the 
mixed phase read 

n a = xh a + (1 - x)n, (47) 
n b = xn h + (1 — x)n, 

with < x < 1. The preferable mixed state for given n a and rib particle densities is the one 
which has the lowest energy 

E Mlx (n a ,n b ) = Min^n {xE N \n a , n b ) + (1 - x)E BCS (n)} , (48) 

where E N , the energy of the normal (unpaired) particles, is given by 



E N (n a ,n b ) 



=<BCS|tf(s = 0)|BCS>= Yl # 

k<P£,a 



(6vr 2 



-5/3 

n a 



-5/3 



(6tt 5 



20tt 2 



'1 — x)n 



x 



5/3 



1 



20tt 2 

rib — (1 — x)n 

x 



+ 



n 



(49) 



5/3 



1 



and E BCS (n) is given by Eq. (|15|). We neglect the interface energy between the normal and 
superconductor components of the mixed phase Fermi gas. For large enough systems the 
surface tension effects are small. Consequently, the mixed phase is energetically favorable 
(and stable) if, and only if, the gain in bulk free energy is larger than the disregarded surface 
energy [3l| . 
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2. On the conditions for the existence of equilibrium between the phases 



The criteria for thermodynamical equilibrium between the components of a mixed phase 
are dictated by the Gibbs conditions 

I. T N = T s , (50) 

II. = 
III. P N = P s . 



These conditions ensure thermal, chemical and mechanical equilibrium, respectively between 
the normal and superconducting phases. Those requirements are enforced by the Maxwell 
construction we employ. Condition I is automatically satisfied since we are considering a 
cold system. As we show below, conditions II and III are reached when we find n and x that 
minimizes the energy. At the minimum we write 

^E^(n a ,n b ) = x^E N (n a ,n b ) + (1 - x)|^E BCS (n) (51) 
= ^P^OE. /£) + A? n a + ffi n h \ + (1 - ^)|^[F BCS (/i) + H = - 

~N , ~N 

Ha + = V, 

where /i = At BCS + /i BCS , fl® = ( - 671 ' ^ 1 and h a = n °'-(i- x ) n _ This resu lt tells us that the sum 
of the chemical potentials between the normal and superconducting phases in equilibrium 
is the same, and not the individual ones. This is because the BCS phase can exchange with 
the normal phase only pairs of a and b species to reach the equilibrium. For the mechanical 
equilibrium, we have 

J^E MIX (n a ,n b ) = E N (n a ,n b ) - E BCS (n) + x^E N (n a , n 6 ) = - (52) 

where we have used the result of Eq. (|51J1 . This is not surprising. We have found already in 
the situation of fixed chemical potentials that coexisting phases have the same thermody- 
namic potential. Since the thermodynamic potential is related to the pressure by F = — P, 
we obtain the desired result P N (/i^,/ib) = P BCS (/i BCS + /i BCS ). Then we have found that 
the mixed phase state satisfies all the necessary conditions to be in thermodynamical equi- 
librium. 



3. The Sarma phase energy 

With Eqs. (J28|) and (|39|). we can write the energy in the Sarma phase as 

E s = F s + ^/i Q n Q (53) 

a 

d 3 k r , „ „ f k2 d 3 k k 2 A 2 



f d k f 



(27r) 3 2m b g 



k>k 2 

Just for completeness, we show in Section IS the expressions for the energy and free energy 
of the Sarma phase and its (expected) results for the limits A — > and A — > A . 
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4- The lowest energy 



Using the conditions for equilibrium between the normal and superfluid phases, Eqs. (15 lj) 
and ()52|) . the mixed phase energy can be rewritten at the minimum as 



E& a ,n 6 ) = x[F N (^,^) +/£ h a + $ h b ] + (1 -x)[F BCS (/i) +H (54) 
= F BCS (/i)+^ n a + $ n b , 

where the n and x dependences in fl^ and pf , are the solutions of Eqs. (J5Tj) and (p)2j) . 
and should be termed and XMm- Both parameters shall depend only, of course, 

m a , m b , n a , n b and a. 

Next we define a function AE(n a , n b ) as the energy difference between the mixed and the 
Sarma phases 

AE(n a , n b ) = E^(n a , n 6 ) - E s (n a , n 6 ) (55) 
= F BCS (/i) +f^n a + n 6 - [F s (/iS, ^) + ^ n a + /if n 6 ] 
= F BCS (/i) - F s (^, ^ + Ua{ ~N _ ^ + nfe( ~N _ ^ 

Now we are able to compare the energy of these two phases for those values of the 
gap supported by the Sarma phase. As we have seen in the situation of fixed chemical 
potentials, the energy difference F S (A ) — F S (A5) is always negative for P F < P FDM (for 
fixed P|). For P b F > P b F 

dmi the mismatch in the Fermi surfaces is large enough so the 
ground state is the normal phase. For the case of a and b fermionic atoms confined in 
an atomic trap we investigate below its ground state in three possibilities that could be 
implemented in current experiments. In the first two cases we vary the interaction through 
the scattering length a and in the last one we fix a and vary the number densities. In most 
of the situations AE(n a ,n b ) has to be evaluated numerically but, as we are going to see 
next, in some limiting cases we can obtain analytical results. 

4-1: Fixed n a and n b with AE(n a , n b , As = 0). When As — > 0, which corresponds to the 
lower limit for the Fermi surfaces mismatch in Eq. (JHHj) . P F — P F 2 = 2y/m a m b A 0j we have 
pS pN ^ gee p^ (|T3|) in the appendix) and consequently [i^ —> Then 

AE(n a , n b , A s = 0) = F BCS (/x) - F N (^, fi) + n a {f a - ft) + n b {f b - tf), (56) 
where = (67r ^f /3 ' . 

For very large asymmetry between the number densities n a and n b , xuin approaches 
1 which implies — > fi^ so the mixed phase would also be reduced to normal phase 
resulting in AE(n a , n b ^> n a , A s = 0) = 0. This happens because with large asymmetry the 
Fermi surfaces are separated apart, making pairing energetically disfavored and the preferred 
ground state is the normal phase. 

We prove analytically below that AE(n a , n b , As = 0) is negative for small asymmetry on 
the number densities of the a and b particle species. In this case the Sarma energy can be 
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written as 



E s {n a ,n b , A s = 0) 



(6 



7T 



2\5/3 



207T 2 



(67T 2 n a 



nT , (n a + <Jn) 5 /3 



1 5/3 



20vr 2 M 



1 + 



5M5n 5M 



+ 7T— — + <W/<) 



3m b n a 9m b \n a J 



(57) 



where 5n = n b — n a is assumed to be small, 5n <C n a . 

An upper bound on AE(n a ,n b , A$ = 0) can be obtained by setting the density of the 
superconducting component of the mixed phase n = n a and minimizing with respect to x. 
Then we have: 



AEK,n 6 ,A s = 0) = -(1 - z)(67r 2 ra a ) 1/3 



(67r 2 n a ) 5 / 3 5n 2 
36n 2 m b n 2 a \x 



MA 2 {P$) 



2tt 2 



(58) 



+ 



1-1 



0(rC 4/3 5n 3 ) 



whose minimization yields 

X %Min 

Thus, the energy difference read 



(6vr 2 n a ) 4 / 3 5n 
lSMm b A 2 (p F )V a 



(59) 



AE(n a ,n 6 ,A s = 0) = -(6n 2 n a ] 



1/3 MA 2 (p F 



27T 2 



XMi 



< 0. 



(60) 



The limit 8n — > implies £Mm ~~ > 0, and E MIX is entirely in the BCS phase. Numerical 
calculations show that the upper bound above is close to the actual minimum. Eq. ()59|) 
shows that the mixed phase is energetically favored compared to the Sarma phase in one 
extreme of the window in Eq. 



4-2: Setting the Sarma gap in its maximum value: AE(n a ,n b , As = A ). This situation 
corresponds to the upper limit for the Fermi surfaces mismatch in Eq. (|38|). P p 2 — P F 2 = 
4y/m a m b A . In this case, As — > A , F s — > F BCS (see appendix), n a = n b = n, + /if) — > 
(/i BCS + fif cs ) and the energy difference in Eq. (|55[) is 

AE(n a = n b = n)= n[$ + ffi - (^ cs + ^ cs )} = 0, (61) 

since for E MIX (n a , n b ), from Eq. (joTj) . we have p% + p% = fj,^ cs + fj^ cs . This means that in 
this situation both the mixed and the Sarma phases are in the BCS state, having the same 
energy. 



4-3: Fixed scattering length varying the number densities with AE(<5n). For intermediate 
values of P F — Pp and for fixed scattering length a (still satisfying the constraint in Eq. 
(|55jl ). the difference AE(n a ,n&) interpolates between the extremes in J^.l and 4-2, as Fig. JSJ) 
exemplifies. We find that for all reasonable values of the parameters (that is, where the mean 
field analysis should apply), and for fixed particle numbers n a and n b , the mixed phase has 
a smaller energy than the Sarma phase. 
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8n/ n 




FIG. 5: Figure on the left: Energy of the Sarma (dashed curve) and mixed phase (solid curve) 
relative to the normal phase as a function of 5n/n where n = (ha + ^b)/2 with Mb = 7 Ma- 
The values of ha and ns were obtained from the Sarma phase for a fixed a, and \a\pA = 0.59, 
0.63 < \a\pB < 0.65. Figure on the right: The fraction of normal state x as a function of bnjn 
with the same masses and scattering length. 



C. Fixed overall density 

We discuss now another case that happens when the overall number of particles N = 
n a + nb is fixed but conversions between the a and b species, a <-> b, are allowed to occur. 
This situation is relevant, for instance, in astro-particle physics. The core of neutron stars 
is believed to be formed by high density quark matter, where different quarks interchange 
their flavors because of weak interactions. In this case, the thermodynamic function to be 
minimized is Li 



K(N = n a + n b ) = E MIX (iV) - xSfx N <5n N - (1 - x)Sfj BCS Sn BCS (62) 
= E MIX -5/i N K-n a ), 



-,N .-.N , T ,,BCS ,,BCS 



where <5/i N = Mb 2 Ma , 5n N = rib — n a , <5/i BCS = — — ^ — > Sn BCS = and we have made use 
ofEq. g7j). 

The particles in the normal component of the mixed phase are accommodated in such a 
way that its energy is given by 



V 7 207T 2 



(j - zfl* (f + zfl* 
m a m b 



(63) 



where N = n a + rib — — 2 ^ x ^ n is the fixed overall number density in the normal fraction 
of the mixed phase. Minimization with respect to z yields 

E»(JV) = (fa ^ )V3 \ (64) 

Then, the most favored mixed state for fixed total density N is the one with the smallest 
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"K energy" : 



K M1X (N = n a + n b ) = Mm x , n { x 



(67T 



2\5/3 



20tT 2 (m 3/2 + m 3/2 )2/3 



N-2{\-x)n 



+(l-x) 



(67r 2 n) 5 / 3 MA2(r2)(67r 2 n) 1 / 3 



20vr 2 M 



2tt 2 



X 



<5/i N (n b - n tt ) ^ . (65) 



5/3 



The two interesting limits that should be evaluated in the situation of fixed overall density 
are m a = m b and m;, ^> m a . 

I: m a = nit,. We set £/i N = to ensure the equilibrium in the conversions between 
the a and b species in the (normal) subsystem of the mixed phase 2 . The chemical 

-2/3 -2/3 

equilibrium implies ^ = Thus, we obtain 

m„ nit 



n a — (1 — x)n + (— ) 3/2 [n b - (1 - z)n]. 



(66) 



Then, since m a = m b this relation gives n a = n b , as it should be. It is found numerically 
that x = and n = n a = n b , favoring the BCS state with energy 



K Mlx (m n 



m b 



(6n 2 n a f/ 3 m a A 2 (n a )(6vr 2 r ia ) 1 /3 



107T 2 m a 



47T 2 



(67) 



II: m b ^> m a . Plugging n a from Eq. (jo^j) in Eq. (jB3j) we find, in the limit m b ^> m a , that 



A' 



MIX 



. 1 (67T 2 ) 5 / 3 

[m b > m a = -^t^T- 2 K - 

(67r 2 n) 5 / 3 



(1 — x)n] 



5/3 



20vr 2 M 



AfA 2 (n)(67r 2 n) 1 / 3 
2^ 



(68) 
(69) 



This equation is minimized when x — 1, giving 



X MIX (m fe » m a ) 



(67r 2 n 6 



,5/3 



207r 2 m b 



(70) 



This is exactly the energy of the normal state with only b particles. It is worth to mention 
that, in a system of fixed overall density with intermediate values of the particles masses 
(different from these two limits we have considered), again the Sarma phase can not be the 
ground state. Since for Sp^ = <5/i BCS = the condition for its existence (Eq. (|38jl) is not 
satisfied. 

Then we have seen that the Sarma phase is always unstable against the mixed phase in 
all situations we have considered. 



2 The same equilibrium is established, of course, in the superfluid component of the mixed phase, i.e. 
«5 M BCS = 0. 
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IV. CONCLUSIONS 



We have studied cold fermionic gases composed by two particle species whose Fermi 
surfaces or densities do not match. We found the ground state of such a systems in three 
physically different situations, namely fixed chemical potentials, fixed number densities and 
fixed overall density. In all cases the extreme limits are the normal and the BCS states, 
depending on the asymmetry in the particle densities or masses. 

For fixed (and different) chemical potentials the system is found to be in the BCS state 
until a point where the asymmetry between the Fermi surfaces is too big so there is a 
first order phase transition to the normal phase. We have investigated the Sarma phase, 
recently proposed as the ground state of asymmetrical fermion systems, employing canonical 
transformation, in mean-field approximation. For fixed chemical potentials, we found that 
the Sarma gap corresponds to an extremum of the free-energy F s . It never represents its 
minimum. 

In the situation where the fixed parameter is the number density, we have found that 
the ground state is a (inhomogeneous) mixed phase in real space, formed by islands of 
asymmetric normal state immersed in a sea of the symmetric BCS state. Again, if the 
particles a and b have the same density, the system is entirely in the BCS state, whereas when 
the asymmetry is above some limiting value, the system is found to be in the normal state. 
We also showed that the two components of the mixed phase state are in thermodynamical 
equilibrium in complete agreement with the Gibbs conditions. 

The third case considered is the one with fixed overall number density, relevant for physics 
of high density quark matter. The two pertinent limits for the masses of the species a and 
b in this case have been considered and we have obtained the favored ground state in these 
situations. 

We have suggested in Ref. (T^J that the mixed phase could be used for detecting su- 
perfluidity in atomic traps. An optical atomic trap should be prepared having an density 
asymmetry in the two species of fermionic atoms. This could be done in current experiments 
|32| , by loading the trap in a way that one of the two hyperfine states of lithium atoms (for 
instance) is more populated than the other. So the islands of normal matter in a sea of BCS 
phase would become apparent and could be imaged. Experiments with different densities 
would have different segregations of the normal phase that could be located, for instance, 
in a corner of the apparatus, leading to a estimative of the "size" of the gap, which is an 
up-to-date open problem. 

Since we do not know how big the surface tension effects in finite systems could be, the 
plans for the future include the study of the magnitude of the interface energy between the 
superfluid and normal phases and also the investigation of the ground state of (asymmetric) 
neutral quark matter. 

We hope that our work will stimulate new experiments with this interesting asymmetrical 
fermion superfluids. 
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V. APPENDIX: EXACT EXPRESSIONS AND LIMITS FOR THE FREE- 
ENERGY AND ENERGY IN THE SARMA PHASE 



The analytical expressions for the free-energy and the energy in the Sarma phase (at the 
minimum), can be written as 



A 5 



F s = 



d 3 k 



q ' Jk<ki (27r) 3 k 

k>k 2 y ' 



ft, p 



30tt 2 M 2tt 2 



Mk F . 9 
-A 2 




(71) 



with e b given by 



k2 d 3 k 



kl (2tt) 



-F b 



\.n 2 m b 



1 3 3 

-(fc 2 — fc-J — (A: 2 — fc-J 



(72) 



The limit A — > results in fci = Pp, hi = P F and Eq. ()71)1 yields the (expected) normal 
free-energy 



F N (/i a ,/i 6 )=F s (A = 0) 
1 



307T 2 



Pf P; 



b5 



(73) 



30tt 2 



(2m a n a ) 5/2 (2m b /i b ) 5/2 



m b 



When the Sarma gap has its maximum value, fci = &2 and the Sarma free-energy reduces to 
the BCS free-energy (Eq. fity) 



^BCS 



(A*) 



ft, pi 



Mfc F 2 
'30tt 2 M 2tt 2 °' 
The analytic expression for the energy in the Sarma phase is 



a 

" k2 k 2 dk k 2 



z.- 5 

h,p 



k 3 A 2 



20tt 2 M 8vr 2 /i 



7T 



\a\kt 



(74) 



(75) 



2tt 2 AM 



&2 L2 



fc 2 dfc fc 2 
:i 2tt 2 2m 6 ' 



Again, with the limit A — > it is very easy to get the normal energy, 



E N (n a ,n b ) = E S (A = 0) 



20vr 2 



&5 



P£ 5 P> 



m a m b 



(67T 



2\5/3 



20tt 2 



n, 



5/3 



5/3 



m b 



(76) 



The BCS energy is recovered in the limit A — > A . In this case the last two terms in Eq. 
(J75J) vanish and we obtain 



E BCS (/x) = E S (A = A ) 



1,5 
i%p 



+ 



P A 2 



20tt 2 M 8vr 2 /i 



7T 



(77) 
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